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Abstract.  We  study  the  detection  of  very  weak  time-periodic  magnetic  signals  via 
a  double-junction  (dc)  Superconducting  Quantum  Interference  Device  (SQUID).  The 
device,  represented  by  two  coupled  nonlinear  differential  equations  for  the  quantum 
mechanical  junction  phases,  admits  of  long-time  static  or  oscillatory  solutions,  the 
transition  being  controllable  by  experimental  parameters.  Signal  detection  is  optimal 
when  the  device  is  “tuned”  to  the  onset  of  the  oscillatory  solutions;  i.e.,  when  the 
minima  in  the  2D  potential  function  disappear.  Modeling  the  device  via  a  derived 
input-output  transfer  characteristic  yields  a  response  (quantified  via  the  signal-to-noise 
ratio  at  the  signal  frequency)  in  good  agreement  with  recent  experiments.  We  also 
present  some  preliminary  results  pertaining  to  coupled  dc  SQUIDs,  taking  into  account 
the  (incoherent)  thermal  noise  sources  in  each  of  the  Josephson  junctions.  As  might 
be  expected  from  previous  work,  coupling  the  SQUIDs  and/or  summing  their  outputs 
enhances  the  response  to  the  external  signal. 


I.  INTRODUCTION 

Superconducting  Quantum  Interference  Devices  (SQUIDs)  [1]  are  the  most  sen¬ 
sitive  detectors  of  magnetic  fields;  however,  they  are  quite  vulnerable  to  environ¬ 
mental  noise,  thermal  noise  in  the  junctions,  or  noise  coupled  in  from  the  biasing 
and  readout  electronics.  Operated  in  the  conventional  flux-locked  mode  wherein 
the  device  is  kept  “locked”  to  an  operating  point  in  the  linear  regime  of  its  trans¬ 
fer  characteristic  via  feedback  electronics,  a  very  small  amount  of  noise  is  often 
sufficient  to  lose  the  operating  point... the  so-called  “slew-rate”  problem.  Allowing 
the  SQUID  to  operate  as  a  free-running  nonlinear  dynamic  device  can  enhance  the 
dynamic  range  if  nonlinear  phenomena  are  carefully  utilized  to  alleviate  noise  prob¬ 
lems.  The  objective  in  this  case  is  not  necessarily  an  enhancement  of  the  SQUID 
response’s  SNR,  rather  it  could  be  an  effective  lowering  of  the  noise  floor  of  the 
readout  system,  resulting  in  an  increase  in  the  entire  system’s  SNR,  signal  detection 
probability,  and  dynamic  range. 

Recent  experiments  and  calculations  have  explored  the  basic  Stochastic  Reso¬ 
nance  (SR)  effect,  and  some  variations,  in  a  single-junction  (rf)  SQUID  [3,4].  This 


20000721  037 

_ HffiKJ  QUALITY  jBSEBeSED-'* - - 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMBNo.  0704-01-0188 


The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources 
gathering  arid  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection 
of  information,  including  suggestions  for  reducing  the  burden  to  Department  of  Defense,  Washington  Headquarters  Services  Directorate  for  Information  Operations  and  Reports 
(0704-0183),  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington  VA  22202*4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be 
subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a  currently  valid  OMB  control  number 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS, _ 


5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 
0601 153N 


5f.  WORK  UNIT  NUMBER 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 
Approved  for  public  release;  distribution  is  unlimited. 


13.  SUPPLEMENTARY  NOTES 


14.  ABSTRACT  This  paper  studies  the  detection  of  very  weak  lime-periodic  magnetic  signals  via  a  double-junction  (dc) 
Superconducting  Quantum  Interference  Device  (SQUID).  The  device,  represented  by  two  coupled  nonlinear  differential  equations 
for  the  quantum  mechanical  junction  phases,  admits  of  long-time  static  or  oscillatory  solutions,  the  transition  being  controllable  by 
experimental  parameters.  Signal  detection  is  optimal  when  the  device  is  "tuned"  to  the  onset  of  the  oscillatory  solutions;  i.e.,  when 
the  minima  in  the  2D  potential  function  disappear.  Modeling  the  device  via  a  derived  input-output  transfer  characteristic  yields  a 
response  (quantified  via  the  signal -to-noisc  ratio  at  the  signal  frequency)  in  good  agreement  with  recent  experiments.  Also  presented 
are  some  preliminary  results  pertaining  to  coupled  dc  SQUIDs,  taking  into  account  the  (incoherent)  thermal  noise  sources  in  each  of 
the  Joscphson  junctions.  As  might  be  expected  from  previous  work,  coupling  the  SQUIDs  and/or  summing  their  outputs  enhances 
the  response  to  the  external  signal. 

Published  in  Proceedings  of  Stochastic  and  Chaotic  Dynamics,  Amblcside,  UK,  16-23  August  1999. 

^SUBJECT  TERMS 
Navigation 

SQUID  noise 

stochastic  resonance  signal  detection 

17.  LIMITATION  OF  Il8.  NUMBER 
ABSTRACT  OF 

PAGES 

UU  13 

« 

Standard  Form  298  (Rev.  8/98) 

Prescribed  by  ANSI  Sid.  Z39.18 


19a.  NAME  OF  RESPONSIBLE  PERSON 
Bulsara,  Adi  R.,  D364 _ 

19B.  TELEPHONE  NUMBER  (Include  area  code) 
(619)  553-1595 


16.  SECURITY  CLASSIFICATION  OF: 

a.  REPORT 

b.  ABSTRACT 

c.  THIS  PAGE 

U 

U 

U 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 
ONR 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Space  and  Naval  Warfare  Systems  Center 

53560  Hull  Street 

San  Diego,  CA  92152-5001 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 
Office  of  Naval  Research 
800  North  Quincy  Street 
Arlington,  VA  22217-5660 


5d.  PROJECT  NUMBER 
MA19 

5e.  TASK  NUMBER 


3.  DATES  COVERED  (From  -To) 
5a.  CONTRACT  NUMBER 


2.  REPORT  TYPE 
Professional  Paper 

4.  TITLE  AND  SUBTITLE  ~  " 

Noise-Mediated  Cooperative  Behavior  and  Signal  Detection  in  dc  SQUIDs 


6.  AUTHORS 
Inchiosa,  Mario  E. 
Bulsara,  Adi  R. 


1.  REPORT  DATE  (DOMM-YYYY) 
08-1999 


device  consists  of  a  single  Josephson  junction  shorted  by  a  superconducting  loop, 
and  coupled  to  the  target  signal  via  cooled  pickup  coils.  The  dynamics  are  in  gen¬ 
eral  multistable,  with  the  magnetic  flux  through  the  superconducting  loop  being 
quantized  in  units  of  the  flux  quantum  4>0  =  h/2e.  In  the  presence  of  the  junction, 
the  magnetic  flux  $  through  the  loop,  in  response  to  an  applied  time-dependent 
magnetic  flux  4>e,  evolves  according  to  the  dynamics  [1], 

.  tlx  =  -U\x).  (1) 

[Note  that  in  this  paper  the  letter  x  (possibly  subscripted)  represents  a  magnetic 
flux  normalized  to  the  flux  quantum;  e.g.,  x  =  d>/d>0.]  In  (1),  x  represents  the 
(normalized)  magnetic  flux  through  the  loop,  the  dot  denotes  time-differentiation, 
and  the  potential  function  U (a:)  =  |(x  —  xe)2  —  cos  2irx  is  multistable  when  /3S  > 
[3 sc  (and  monostable  otherwise);  fisc  —  1  for  xe  —  0.  xe  —  x,-  +  xo  is  the  externally 
applied  flux,  xo  is  a  dc  magnetic  flux,  and  x,-  —  r/(t)  +  y(t )  is  an  ac  magnetic 
flux  consisting  of  a  deterministic  component  rj{t)  and  a  background,  non-thermal 
noise  y(t)  (fed  in  from  the  bias  or  readout  electronics  or  the  environment).  The 
noise  is  usually  modeled  as  Gaussian  and  exponentially  correlated,  with  correlation 
time  much  less  that  the  SQUID  time  constant  tl,,  which  is  typically  the  smallest 
time-scale  of  the  problem  (here  we  are  assuming  the  SQUID’s  thermal  noise  to  be 
negligible).  The  inertial  (capacitative)  term  in  the  dynamics  is  absent  because  in 
most  applications  one  incorporates  a  shunt  resistance  across  the  junction  [1];  this 
removes  hysteresis  in  its  voltage-current  characteristic. 

The  signal  detection  properties  of  the  rf  SQUID  have  been  extensively  analyzed 
for  the  case  of  a  multistable  potential  function.  In  the  presence  of  a  known  time- 
periodic  bias  signal,  an  unknown  dc  “target”  magnetic  flux  Xo  renders  the  potential 
asymmetric,  and  all  harmonics  of  the  (known)  sinusoidal  bias  flux  y(t)  appear.  An 
SR.  effect  occurs  at  every  harmonic  and  the  spectral  amplitudes  are  dependent  on 
the  magnitude  of  the  target  dc  signal.  For  a  potential  whose  minima  all  come  in 
symmetric  pairs  (e.g.  a  symmetric  double-well  potential),  only  the  odd  multiples 
of  the  signal  frequency  appear  in  the  power  spectral  density  (PSD)  of  the  motion. 
Hence,  the  appearance  of  the  even  multiples  of  the  bias  signal,  together  with  the 
change  in  their  spectral  amplitudes  can  be  used  as  a  means  to  detect  and  quan¬ 
tify  the  target  signal  [4].  In  complex  operating  scenarios,  if  a  “clean”  window  is 
known  to  exist  in  a  particular  regime  of  the  frequency  spectrum,  we  can  choose 
the  frequency  of  ?/(f)  to  carry  out  our  detection  in  that  window.  This  technique 
makes  no  attempt  to  eliminate  low  frequency  noise  or  to  enhance  the  output  SNR 
in  the  conventional  sense;  it  merely  shifts  the  detection  to  a  more  acceptable  part 
of  the  frequency  spectrum.  This  idea  affords  a  novel  technique  that  can  be  applied 
to  SQUIDs  (particularly  high  temperature  SQUIDs)  as  well  as  other  nonlinear  dy¬ 
namic  sensors  that  arc  noise-constrained  at  certain  frequencies. 

In  the  monostable  regime,  it  is  convenient  to  characterize  the  SQUID  [6]  as  a 
memoryless  static  nonlinearity,  a  /ion-dynamical  system  characterized  only  via  an 
input-output  transfer  characteristic  (TO).  This  characteristic  is  found  by  setting 


the  left  hand  side  of  (1)  to  zero,  solving  the  resultant  transcendental  equation  for 
x{t)  and  plotting  the  “shielding  flux”  xs(t)  =  x(t)  —  xe(t)  vs.  the  input  flux  #,(£). 
Having  obtained  the  TC,  one  can  analytically  apply  the  sinusoidal  input  signal  plus 
Gaussian  noise  to  it,  and  compute  the  autocorrelation  function  {xs(t)xs(t  +  h )} 
and  its  associated  power  spectral  density  [6],  whence  we  compute  the  SNR  at 
the  fundamental  of  the  applied  signal.  On  a  3D  plot  of  SNR  vs.  dc  bias  x0  and 
nonlinearity  j3s,  one  observe^  distinct  regimes  of  SNR  crests  and  troughs.  This  form 
of  SNR  maximization  differs  from  SR  insofar  as  one  obtains  the  optimal  response 
by  tuning  control  parameters  (in  tills  case,  the  dc  bias  flux  and  the  nonlinearity 
parameter),  rather  than  the  input  noise  strength.  The  theoretical  results  are  well 
supported  by  experiments,  and  the  response  of  such  static  nonlinearities  to  more 
general  (e.g.  wideband)  signals  has  also  been  investigated  in  a  recent  publication 
[7]  using  information-theoretic  distance  measures  to  characterize  the  output. 

With  this  preamble,  we  now  turn  to  our  recent  work  on  the  two-junction  or  dc 
SQUID  [8,9].  Unlike  the  rf  SQUID  described  above,  the  dc  SQUID  is  characterized 
by  a  2D  potential  function  and  can  exhibit  spontaneous  oscillatory  solutions  under 
the  appropriate  system  preparation.  We  first  describe  the  single  dc  SQUID,  and 
then  outline  some  ongoing  research  on  the  response  of  two  flux-coupled  dc  SQUIDs. 

II  NONLINEAR  AMPLIFICATION  IN  A  DC  SQUID 

The  dc  SQUID  consists  of  two  Josephson  junctions  inserted  into  a  superconduct¬ 
ing  loop  [1];  we  assume,  for  convenience,  that  the  insertion  is  symmetric.  Con¬ 
ventionally,  the  voltage  measured  across  the  Josephson  junctions  is  taken  as  the 
SQUIDs  “output.”  Instead^  we  take  the  circulating  current  Is  (experimentally  mea¬ 
sured  via  the  associated  “shielding  flux”)  as  output.  This  setup  was  developed  for 
studying  SR  in  dc  SQUIDs  operating  in  a  hysteretic  regime.  However,  much  higher 
output  signal  strengths  and  SNRs  were  discovered  by  using  dc  bias  currents  large 
enough  to  take  the  device  beyond  the  hysteretic  regime  into  the  regime  of  oscilla¬ 
tory  solutions.  These  higher  input-output  gains  result  from  the  rapid  change  of  Is 
with  <I>C  close  to  where  the  dynamics  change  from  static  to  oscillatory. 

In  the  presence  of  an  external  magnetic  flux  <f>e,  one  obtains  [1]  a  loop  flux 
consisting  of  the  (geometrical)  component  d>c.  together  with  a  contribution  arising 
from  the  induced  circulating  or  shielding  current  Is  that  tends  to  screen  the  applied 
flux: 


<I>  =  (I>e  +  LIS,  (2) 

L  being  the  loop  inductance.  The  Josephson  currents  in  each  arm  of  the  “interfer¬ 
ometer”  are  /q  sin  and  /0  sin  J2,  with  the  junctions  assumed  to  be  identical  with 
critical  currents  /o,  and  with  <5i>2  being  the  quantum  phases.  The  wave-function 
must  remain  single-valued  around  the  SQUID  loop,  leading  to  the  phase  continuity 
condition, 


$2  —  <5i  —  2nn  —  2ir$/%y  (3) 

n  being  ail  integer,  and  4>o  =  k/2c  the  flux  quantum.  Combining  (2)  and  (3)  and 
setting  n  —  0,  we  find  for  the  circulating  current  Is, 

0~r  =  Si  —  82 
•  /o  $0 

where  0  =  2nLIo/$o  is  the  nonlinearity  parameter.  In  the  absence  of  noise  and 
the  target  magnetic  flux  (taken  to  be  sinusoidal  in  this  work),  we  can  use  the  RSJ 
model  to  write  down  equations  for  the  currents  in  the  two  arms  of  the  SQUID  via  a 
lumped  circuit  representation  [1];  expressed  via  the  Josephson  relations  Si  —  2 eVi/h 
linking  the  voltage  and  the  quantum  phase  difference  across  the  junction  i,  these 
equations  take  the  form, 


tSl  —  ^  -  Is  —  /0sin£i,  tS-2  =  ^  +  Is  —  losmS-z, 


(5) 


where  r  =  tific/2,  li  being  the  normal  state  resistance  of  the  junctions.  The  dc  bias 
current  h_,  is  applied  symmetrically  to  the  loop.  I11  experiments  [8],  the  bias  current 
and  applied  flux  are  externally  controllable.  This  is  a  critical  point,  since,  as  will 
become  evident  below,  it  permits  us  to  manipulate  the  shape  of  the  2D  potential 
function  that  characterizes  the  SQUID  dynamics  and  thereby  the  input-output  TC 
that  governs  the  response  (note  that  the  solutions  of  (5)  can  be  oscillatory,  even 
in  the  absence  of  external  inputs).  Rescaling  the  time  by  r//o,  one  can  write  the 
above  in  the  form  Si  —  —jf.  with  the  2D  potential  function  defined  as 


U{SUS2)  -  -  cos  Jx  -  cos  S2  -  J(Si  +  S2)  +  (20)  ^  -  S2  -  2n$ex)2,  (6) 

where  we  introduce  the  dimensionless  bias  current  J  =  /(,/(2/o)  and  normalized 
applied  flux  (I>er  =  d>c./4>o,  taken  to  be  solely  dc  for  the  moment. 

The  SQUID’s  Josephson  junctions  are  edways  in  a  superconducting  state  when 
the  potential  (6)  has  stable  minima;  it  is  readily  apparent  that  the  symmetry  of  the 
potential  and  the  depth  of  the  minima  are  controlled  by  the  adjustable  parameters 
J  and  (I>cr.  This  configuration  (including  the  problem  of  thermal  activation  out  of 
the  stable  states  of  the  potential)  lias  been  discussed  in  the  literature  [10].  After 
a  brief  transient,  the  phase  angles  A12  achieve  constant  steady-state  values  and 
one  obtains  the  conditions  for  the  minima  via  Si<2  =  0.  This  leads  to  the  current 
equations 


=  /0(sinJ’i  +  sin  £2),  2IS  —  I0(smS2  —  sin  Si).  (7) 

01  course,  these  equations  may  <dso  be  written  down  by  applying  Kirchoff’s  laws 
directly  to  the  lumped  circuit  representation  of  the  SQUID.  Using  the  phase  conti¬ 
nuity  relation,  we  are  finally  able  to  write  down  a  transcendental  equation  for  the 
circulating  current  Is: 


(8) 


k 

/o 


sin  fyr^e!-  +  jjA  cos 


Sill 


-1 


+  n$ex  +  ^ry- 
LIq 


Equation  (8)  may  be  solved  numerically  for  the  circulating  current;  the  ensuing  TC 
is  periodic  in  the  applied  flux  <3>ex.  and  possibly  hysteretic,  with  the  hysteresis  loop 
width  controlled  by  the  bias  current  J.  For  J  —  0  one  obtains  hysteresis  for  any 
nonlinearity  fi;  for  0  <  J  <  1,  hysteresis  occurs  over  some  range  of  /?.  It  is  most 
important  to  note  that  equations  (7)  and,  therefore,  (8)  are  valid  only  when  the 
potential  has  stable  minima.  In  this  regime,  the  externally  applied  bias  current  is 
matched  by  the  sum  of  the  junction  supercurrents.  When  this  balance  is  exceeded, 
a  finite  voltage  V  (corresponding  to  a  normal  loop  current  V/R)  appears  across  the 
device.  The  maximum  applied  current  (before  the  device  enters  its  “voltage  state”) 
is  clearly  2/0  in  the  absence  of  an  applied  external  flux.  However,  in  the  presence 
of  the  external  flux,  one  must  compute  the  the  critical  applied  current  at  which 
the  voltage  state  appears.  In  this  regime,  one  obtains  oscillatory  solutions  for  the 
phases  Si  (modulo  27t)  and  the  circulating  current  Is,  reminiscent  of  self-excited  or 
relaxation  oscillators  that  are  encountered  in  systems  with  negative  damping  [11]. 
For  the  system  at  hand,  this  behavior  may  be  traced  to  the  description  in  terms 
of  the  two  coupled  first-order  differential  equations  (5),  together  with  a  non-zero 
applied  current  J .  In  fact,  it  may  be  shown  that,  only  in  the  superconducting 
regime  (where  the  potential  has  stable  minima),  the  SQUID  dynamics  (5)  may  be 
reduced  to  the  ID  form,  in  terms  of  the  normalized  flux  variable  x: 

rsx  —  —x  —  xc  —  —  sin  ttx  cos  Z,  (9) 

where  we  set  rs  =  ^  and  Z  =  irx  +  arc-sin  [ J  +  ^f{x  —  xe)j.  Equating  the  right 
hand  side  of  (9)  to  zero  and  solving  (numerically)  for  the  x  vs.  xe  TC  yields  a  curve 
identical  to  that  obtained  from  the  transcendental  form  (8),  after  we  express  Is  in 
terms  of  x.  Note  Unit,  in  this  regime,  the  potential  (6)  can  readily  be  transformed 
into  a  single- variable  potential  U(x),  whose  gradient  yields  (up  to  a  multiplicative 
constant)  the  negative  of  the  right  hand  side  of  (9);  this,  of  course,  is  to  be  expected. 

We  are  interested,  primarily,  in  operating  the  SQUID  in  the  regime  that  yields 
the  optimal  response  to  a  target  input  signal;  further,  we  would  like  to  (as  outlined 
in  the  preceding  section)  obtain  this  response  only  by  adjusting  the  system  parame¬ 
ters,  without  hciving  to  alter  the  input  SNR.  From  a  computational  standpoint,  it  is 
certainly  more  convenient  to  work  in  the  non-hysteretic  regime  wherein  the  output 
variable  (in  this  case,  the  circulating  current  Is)  is  single-valued  in  the  input.  Then, 
one  simply  constructs  the  TC!  and  applies  the  signal  and  noise  to  the  input,  as  has 
been  done  in  our  earlier  work  [6].  In  the  hysteretic  regime,  on  the  other  hand,  one 
must  integrate  the  coupled  equations  (5)  and  compute  power  spectra  by  averaging 
time  series;  this  can  prove  cumbersome  in  some  cases,  and  is  certainly  dependent  on 
computing  power.  For  the  rf  SQUID,  outlined  in  the  preceding  section,  hysteresis 
is  removed  by  decreasing  the  nonlinearity  parameter  below  a  critical  value.  In  the 


FIGURE  1.  (a)  Transfer  characteristic:  time-averaged  circulating  current  (see  text)  I,/Io  vs. 
applied  Uux  4>ex  for  fl  —  0.4.  (b)  Contour  plot  showing  theoretically  predicted  output  SNR  (taken 
at  target  signal  frequency),  in  the  oscillatory  solutions  regime,  vs.  bias  parameters  ^bias  and  J. 
SNR  scale  (black-to-white)  corresponds  to  -28  dB  to  32  dB,  with  contour  lines  spaced  5  dB  apart 
(SNR  values  of  -28  dB  or  less  are  represented  by  black;  a  “speckled”  appearance  in  the  black 
regions  of  the  /?  —  0.4  plot  is  due  to  limited  numerical  precision). 


current  situation,  however,  one  may  obtain  hysteresis  for  any  /?,  with  the  applied 
current  J  controlling  the  regime  of  hysteresis.  In  experiments  and  simulations,  one 
obtains  the  best  response  when  the  TO  is  non-hysteretic;  this  might  be  expected 
since,  in  general,  the  optimal  response  occurs  in  the  linear  or  quasi-linear  regime. 
The  bulk  of  our  work  is,  therefore,  carried  out  in  the  non-hysteretic  regime,  with 
fi  ~  1  used  in  most  of  our  simulations.  In  this  regime,  for  a  range  of  applied  fluxes 
around  4>ea.  =  0.5  the  SQUID  admits  of  oscillatory  solutions  and  the  potential 
ceases  to  have  minima,  with  the  conservation  relations  (7)  violated. 

We  now  examine  the  dc.  SQUID  in  the  regime  of  oscillatory  solutions  in  some 
greater  detail.  These  solutions  arise  when  the  bias  current  J  is  increased  to  the 
point  where  the  first  of  the  conservation  relations  (7)  cannot  hold  without  the 
presence  of  an  additional  (normal)  current  term.  Let  $ea.i  be  the  critical  flux  at 
which,  as  (l>ei-  is  increased  from  0  to  1/2,  the  potential’s  minima  disappear  and 
oscillatory  solutions  replace  static  ones.  For  increasing  from  1/2  to  1,  the 
oscillatory  solutions  disappear  at  <i>ca.2  =  1  —  dhx-i-  Recent  calculations  [9]  predict 
the  critical  external  dc  flux  <l>c.rl ,  for  a  given  bias  current  J .  The  oscillation  period 
tends  to  infinity  very  close  to  4>exi,  decreases  as  d>ear  — >■  1/2,  and  increases  again  as 
<1>„.  —>  <I>er2.  The  interval  dhu-i)  of  the  oscillatory  solutions  increases  with  J. 

Other  researchers  [12]  have  investigated  the  oscillatory  regime  via  a  computation 


of  the  voltage  across  the  Josephson  junctions;  the  Josephson  relations  relate  this 
voltage  (which  can  be  experimentally  measured)  to  the  frequency  of  the  oscillations. 
Here,  however,  we  consider  the  screening  current  Is  as  the  variable  of  interest.  Our 
experiments  [8]  are  able  to* measure  the  current,  with  one  important  caveat:  due 
to  the  extremely  high  oscillation  frequency  [predicated  by  the  constant  r  in  (5)], 
the  experiments  actually  measure  the  time-averaged  current  Is.  The  TC  is  then 
simply  a  plot  of  this  /,  vs.  the  applied  flux  $er.  An  analytic  computation  of  Is  has 
not  yet  been  accomplished;  however,  one  may  numerically  solve  the  dynamics  (5), 
compute  the  (oscillating)  Is  and  then  compute  the  time-average.  A  family  of  TCs 
obtained  in  this  manner  is  shown  in  Fig.  1(a). 

From  the  TCs,  we  can  calculate  the  output  SNR  in  response  to  an  applied  flux 
consisting  of  a  time-sinusoidal  signal,  colored  noise,  and  a  dc  component  $bias 
(normalized  to  4>o).  The  input  signal  power  Sm  and  noise  power  N-m  are  taken 
arbitrarily  small  so  that  they  are  transformed  in  a  nearly  linear  manner  by  the 
device,  simply  being  multiplied  by  the  slope  or  gain  G  of  the  TC  at  the  bias  value 
d’bias.  Note  that  the  strongly  nonlinear  response  regime  has  been  studied  in  [6]. 
We  assume  that  the  signal  and  noise  lie  well  within  the  SQUID  bandwidth  t^1.  In 
the  experiment  a  second  SQUID,  operating  as  a  conventional  (flux-locked)  device, 
reads  the  output  of  the  “sensing”  SQUID;  hence  it  introduces  a  noise  floor  Nj 
(taken  equal  to  N-m)  into  the  measurement  of  the  sensing  SQUID.  Taking  into 
account  the  power  coupling  efficiency  e  between  the  SQUIDs,  we  may  write  down 
for  the  output  SNR, 

'  S-mG2c  =  RmG2c 
N-mG2t  +  Nj  Gh  +  1  ’ 

in  Fig.  1(b)  we  show  the  output  SNR  vs.  the  control  parameters  J  and  $bias-  The 
4’bias  values  shown  cover  two  periods  of  the  TC,  and  the  SNR  contour  plots  reflect 
this  periodicity.  The  J  values  shown  cover  the  range  over  which  the  dc  SQUID 
has  both  static  and  oscillatory  states  (depending  on  the  value  of  $t>ias)-  This  range 
is  characterized  by  a  split  in  the  maximal  SNR;  one  observes  “bifurcating  crests” 
separated,  vertically,  by  (I>Cr2— The  SNR  depends  on  the  slope  of  the  TC,  being 
maximal  on  the  segments  with  maximum  slope.  The  troughs  running  just  outside 
the  crests  correspond  to  zero-slope  points  at  minima  and  maxima  of  the  TC.  With 
increasing  fi,  the  range  of  J  values  yielding  both  static  and  oscillatory  solutions 
shrinks.  Also,  for  large  /Fs  the  crests  become  less  evident,  but  the  troughs  remain. 
For  J  >  i,  the  Josephson  junctions  are  normal  and  the  transfer  characteristic  has 
a-  much  reduced  range,  resulting  in  a  minimal  output  SNR. 

The  theoretical  predictions  have  been  shown  [9]  to  be  in  very  good  agreement 
with  experimental  results  [8].  It  is  particularly  gratifying  to  be  able  to  “tune” 
the  response  (i.e.  adjust  the  TC)  to  the  optimal,  by  adjusting  the  experimentally 
accessible  parameters  J  and,  $  bias?  since  the  nonlinearity  /?  cannot  easily  be  adjusted 
after  fabrication.  The  agreement  with  experimental  results  could  be  substantively 
improved  by  a  systematic  consideration  of  coupling  effects  between  the  sensing 


and  readout  SQUIDs.  It  is  also  important  to  study  coupling  effects  as  a  potential 
means  of  reducing  the  effective  noise  floor  of  the  readout;  our  earlier  work  [13]  on 
coupled  systems  indicates  that  substantial  enhancements  of  the  output  SNR  should 
accrue,  with  carefully  applied  coupling.  In  the  following  section,  we  outline  some 
preliminary  results  for  coupled  tie  SQUIDs. 


Ill  COUPLED  DC  SQUIDS 


Motivated  by  the  plethora  of  recent  results  [13]  showing  that  coupling  and  noise 
can,  cooperatively,  enhance  the  response  in  nonlinear  dynamic  devices,  we  consider 
and  present  some  preliminary  results  involving  coupled  dc  SQUIDs.  In  general,  the 
coupling  can  be  globed  or  local,  linear  or  nonlinear.  Here  we  consider  n  SQUIDs 
with  global  linear  coupling.  The  klh  SQUID  (k  —  l,...,n)  is  described  by  the 
(coupled)  dynamics  of  the  phase  angles  5^  (j  —  1,2): 

=  Jt  +  -  sin*,-,  (11) 

Rk  h  k 

where  R  represents  the  circulating  current  (we  have  dropped  the  subscript  s  for 
convenience),  and  J/;  —  R\R (2/0/,-)  the  (normalized)  externally  applied  bias  current; 
I  ole  is  the  critical  current  of  the  junctions,  and  ta  =  HR^-e/  2,  with  Rk  being  the 
normal  state  resistance  of  t[ie  junctions  (we  assume  that  each  SQUID  is  comprised 
of  two  identical  junctions,  with  the  differences  in  /Us  arising  via  different  loop 
inductances  La  ).  The  circulating  current  is  inductively  coupled  (we  assume  equal 
mutual  inductance  coupling  of  strength  M)  to  the  loop  currents  in  the  remaining 
SQUIDs: 
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where  4>tAo  is  the  external  flux  at  the  klh  SQUID.  Noting  that  the  bias  currents  Im 
appearing  on  the  right  hand  side  of  (12)  are,  in  turn,  coupled  to  every  other  bias 
current  (including  R),  we  see  that  we  will  get  an  infinite  nested  series  on  the  right 
hand  side  of  (12)  and,  hence,  in  (11).  For  the  purposes  of  this  paper,  we  consider 
some  simplified  situations,  namely  a  pair  of  (possibly  non-identical)  SQUIDs  and 
n  identical  SQUIDs. 

First  we  consider  the  scenario  of  two  inductively  coupled  dc  SQUIDs  (for  this 
case  global  and  local  coupling  are  of  course  equivalent).  The  circulating  currents 
in  the  two  loops  evolve  via  dynamics  analogous  to  (4): 
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From  (11),  we  may  write  down  equations  for  the  coupled  phase  angles  for  the 
SQUIDs,  after  some  manipulation: 
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where  K12  =  i<-\J hi/ Z/2,  k-21  =  kJlJLI  and  /?,-  =  /?,(1  —  k2),  with  «  =  M/y/L iL2. 
Note  that  /c  can  range  from  zero  to  one. 

Before  introducing  signal  and  noise  terms,  we  note  a  special  case.  If  the  SQUIDs 
are  identical  and  identically  biased  (Ji  =  =  J,  $eio  =  $e2o  =  $eo)>  the  dynamics 

above  reduce  to  the  much  simpler  single  SQUID  equations  of  motion  (5)  with  a 
new  effective  nonlinearity  parameter  f)  =  /i(l  +  k).  Rescaling  time  by  gives 
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which  may  be  treated  as  described  in  the  preceding  section  and  our  recent  work 
[9].  in  the  absence  of  noise,  the  loop  fluxes  evolve  identically,  with  the  system  (15) 
describing  the  evolution  of  each  SQUID.  The  question  of  phase  locking  in  a  similar 
multi-junction,  two  loop  system  was  studied  in  [14].  Note  that  the  general  case  of 
n  globally  coupled,  identical,  and  identically  biased  SQUIDs  may  be  characterized 
rather  elegantly  after  iterative  calculations  starting  with  (11);  for  this  case,  we  are 
able  to  write  the  dynamics,  in  the  absence  of  noise  terms,  in  closed  form  analogous 
to  (14),  whence  we  find  that  the  dynamics  collapses  onto  the  description  of  a  single 
SQUID  with  a  “dressed”  (i  as  in  (15): 
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tc  —  M/L,  from  which  we  recover  (15)  for  n  —  2. 

We  now  consider  the  interesting  situation  of  n  identical  and  identically  biased 
SQUIDs  subject  to  different  (internal)  noise  sources.  The  equations  for  the  phase 
angles  may  be  obtained  directly  from  (11)  by  adding  a  (thermal)  noise  term,  as¬ 
sumed  to  be  Gaussian  white  noise  Fkj(t)  (k  =  1, . . .  ,n;  j  =  1,2),  to  each  equation 
to  model  thermal  noise  arising  in  the  Josephson  junctions.  The  noise  sources  are 
assumed  uiicorrclatcd,  and  each  source  has  autocorrelation  (Fkj(t)Fkj(t  +  s)}  — 
Di;j$(s).  Note  that  even  with  identical  and  identically  biased  SQUIDs,  we  cannot 
collapse  the  dynamics  onto  two  equations  as  in  (16)  because  of  the  noise  terms.  In 
the  absence  of  noise,  however,  the  description  (16)  may  be  used  to  determine  the 


critical  bias  conditions  for  the  onset  of  the  oscillatory  solutions,  as  in  the  preceding 
section. 

We  can  compare  severed  interesting  cases:  a  single  SQUID  (n  =  1),  an  array 
of  uncoupled  SQUIDs  (?i  >  1,  k  =  0),  and  an  array  of  coupled  SQUIDs  (n  >  1, 
k  >  0).  Note  that  if  we  take  cis  the  array’s  output  the  circulating  current  in  one 
member  of  the  array,  for  k  —  0  we  will  clearly  obtain  the  same  results  for  n  =  1 
and  n  >  1;  however,  if  we  take  as  the  output  the  sum  of  the  circulating  currents  of 
all  of  the  array  elements,  these  two  cases  become  distinct. 

We  simulated  the  above  cases  using  (3  =  2  and  an  externally  applied  flux  consist¬ 
ing  of  a  dc.  bias  and  a  sinusoidal  target  signal:  (i\.0  =  4>er  +  Asin(te£)  (k  —  1, . . .  ,n). 
We  bias  the  single  SQUID  such  that  the  TO  looks  approximately  like  that  in 
Fig.  1(a),  panel  J  =  0.1,  having  a  steep  region  around  4>er  —  0.5.  To  effect 
this  biasing  for  our  (3  —  2  SQUID,  we  set  the  bias  current  to  J  =  0.40731,  giving 
a  critical  flux  for  the  onset  of  running  states  of  $ea,i  =  0.495  and  producing  a  very 
steep  region  in  the  TC  from  about  (I>£.ri  —  0.495  to  about  <ber2  —  0.505.  The  gain 
(slope  of  the  TC!)  at  <l>ex  =  0.5  in  the  absence  of  noise  is  approximately  53.1.  We 
then  set  the  bias  current  J  of  the  coupled  array  such  that  the  predicted  gain  (in  the 
absence  of  noise)  equals  the  gain  of  the  single  SQUID.  For  example,  for  k  =  0.1, 
the  bias  current  needed  to  yield  a  gain  of  53.1  in  a  four  SQUID  array  is  J  =  0.4777. 
Note  that  (in  the  absence  of  noise)  we  can  reduce  the  problem  of  predicting  the 
gain  of  the  array  to  one  of  predicting  the  gain  of  a  single  SQUID,  using  the  fact 
that  we  can  reduce  the  n  coupled  SQUIDs  to  a  single  one  with  a  “dressed”  (3. 

The  circulating  current’s  oscillation  frequency  is  an  important  time  scale  to  con¬ 
sider  when  performing  simulations.  The  step  size  of  the  numerical  integration 
should  be  much  smaller  than  the  period  of  this  oscillation.  As  mentioned  in  the 
previous  section,  the  circulating  current’s  oscillation  frequency  is  typically  very 
high,  and  the  electronics  is  only  able  to  read  out  the  time-averaged  value  of  the 
circulating  current.  In  the  simulations,  we  effectively  perform  this  averaging  by 
looking  only  at  frequencies  in  the  power  spectrum  close  to  the  target  signal,  which 
we  set  at  a  value  much  less  that  the  circulating  current’s  oscillation  frequency. 

Since  we  have  chosen  the  bias  currents  to  yield  equal  gains  in  the  single  SQUID 
and  the  array,  we  expect  to  get  the  same  signal  power  out  of  both,  at  least  in  the 
limit  of  weak  signal  and  noise.  Concerning  the  noise,  note  that  each  SQUID  in  the 
array  is  subject  to  the  same  input  signal,  but  each  Josephson  junction  gives  rise  to 
its  own  thermal  noise,  uncorrelated  with  the  other  noise  sources.  Thus,  coupling 
the  SQUIDs  or  summing  their  outputs  will  give  rise  to  an  “averaging”  over  these 
un correlated  noise  sources.  This  should  result  in  an  increased  SNR. 

Let’s  consider  the  ease  of  summing  in  more  detail.  All  n  of  the  SQUIDs  are 
subjected  to  the  same  target  signal,  so  their  responses  at  the  target  signal  will  all 
be  correlated  with  each  other.  If  they  are  perfectly  correlated,  the  amplitude  of  the 
summed  output  will  be  n  times  greater  than  a  single  output,  and  the  summed  signal 
power  will  be  n2  times  greater.  In  contrast,  the  noise  generated  in  the  SQUIDs  will 
be  uncorrelated  from  SQUID  to  SQUID.  The  powers  (rather  than  the  amplitudes) 
ol  random,  uncorrelated  noises  add,  so  the  summed  noise  power  will  be  only  n 


times  greater  than  that  at  a  single  output.  Hence,  in  the  ideal  case  of  perfectly 
correlated  signals  and  uncorrclated  noises,  summing  boosts  the  SNR  by  a  factor  of 
n2 /'a  —  n  compared  to  the  SNR  of  an  individual  output. 

In  practice,  the  summed  output  could  be  measured  in  a  variety  of  ways.  We 
could  individually  measure  the  shielding  flux  of  each  SQUID,  or  we  could  use  a 
network  of  pickup  coils  connected  in  parallel  to  sense  and  sum  the  shielding  fluxes 
of  all  of  the  SQUIDs. 

Certainly  the  network  of  pickup  coils  would  be  much  simpler  to  implement  than 
measuring  each  SQUID  individually.  A  further  advantage  of  measuring  the  summed 
output  (rather  than  measuring  the  individual  outputs  and  then  summing  the  mear 
surements)  would  be  that  summing  would  also  boost  the  signal  power  into  the 
measurement  device,  and  if  the  noise  floor  of  the  measurement  device  is  the  limit¬ 
ing  factor  in  output  SNR,  this  will  increase  the  output  SNR. 

Of  course,  the  summing  network  will  also  unavoidably  cause  a  coupling  between 
the  SQUIDs.  In  fact,  this  sort  of  coupling  is  identical  to  the  global  coupling  case 
we  have  been  considering.  Each  SQUID  feels  a  flux  due  to  all  the  other  SQUIDs 
(but  not  from  itself,  of  course).  But  how  would  this  coupling  affect  the  SNR 
improvement  produced  by  summing?  By  causing  more  correlation  between  the 
SQUIDs,  strong  coupling  could  conceivably  kill  the  improvement.  Our  simulations 
address  this  question  for  an  array  of  four  SQUIDs. 

By  summing  the  circulating  currents  of  4  uncoupled  SQUIDs,  we  get  a  6  dB 
improvement  over  a  single  SQUID.  This  is  what  one  would  predict  from  theory: 
by  summing  4  uncoupled  SQUIDs  one  gets  42  =  16  times  the  signal  power  (12 
dB  increase)  and  just  4  times  the  noise  power  (6  clB  increase)  of  a  single  SQUID, 
yielding  4  times  the  SNR  (6  dB  improvement).  Figure  2  illustrates  this  result  for  a 
target  signal  amplitude  of  A  —  0.000125893  and  frequency  u>  —  7r / 1000  [using  the 
rescaled  time  units  of  equation  (16)].  In  the  output  SNR  plot,  the  uncoupled  array 
(k  =  0)  “summed  output”  scenario  consistently  exceeds  the  output  SNR  of  a  single 
uncoupled  SQUID  by  approximately  6  dB,  over  the  range  of  Josephson  junction 
noise  strengths  plotted.  The  output  signal  power  plot  shows  the  uncoupled  array 
“summed  output”  scenario  exceeding  the  output  signal  power  of  a  single  uncoupled 
SQUID  by  approximately  12  clB. 

When  we  simulate  a  system  with  moderate  coupling  (k  —  0.1)  between  elements, 
we  find  nearly  the  SNR  obtained  by  summing  the  outputs  of  uncoupled  elements. 
Referring  again  to  Fig.  2,  we  find  an  approximately  4.6  dB  SNR  improvement 
with  summing  and  coupling,  while  coupling  boosts  the  SNR  of  an  individual  array 
element  by  about  2.3  dB.  For  low  noise  strengths,  the  coupled  array  summed  output 
signal  power  matches  that  of  the  uncoupled  array,  but  for  higher  noise  the  effect  of 
coupling  is  to  reduce  the  power  somewhat. 

in  this  system,  summing  and  coupling  does  not  seem  to  beat  summing  alone. 
However,  the  important  observation  is  that  moderate  coupling  does  not  seriously 
hurt  the  summing  effect.  This  means  that  in  practice  the  simple  summing  network 
described  above  can  be  used,  and  the  side  effect  of  coupling  will  not  kill  the  SNR- 
improving  effect  of  summing. 


FIGURE  2.  Output  SNR.  ft  and  output  signal  power  S.  Solid  line:  k  —  0,  single  output. 
Dot-dashed  line:  k  =  0.1,  single  output.  Dashed  line:  k  =  0.1,  summed  output.  Dotted  line: 
k  —  0,  summed  output.  All  scales  are  in  decibels,  i.e.  101og10(-). 


The  results  shown  in  Fig.  2  are  for  a  “weak”  signal,  meaning  that  A  <  ($ex2  — 
<l>erl  )/2.  Thus,  the  response  to  the  target  signal  is  nearly  linear  provided  the  noise  is 
also  weak.  We  have  obtained  qualitatively  similar  results  for  strong  signals  having 
A  (Ihr 2  —  (l>exi)/2,  but  in  this  case  the  response  nonlinearity  causes  the  output 
signal  power  (summed  or  not)  to  be  slightly  higher  with  k  —  0.1  than  with  k  —  0, 
provided  the  noise  is  weak.  As  the  noise  level  is  increased,  the  fact  mentioned 
above  concerning  a  slight  cpupling-related  power  reduction  at  higher  noises  causes 
the  curves  for  h  =  0.1  and  k  =  0  to  cross,  with  the  coupled  case  falling  below  the 
uncoupled  case. 

In  simulations  with  ic  =  0.5  (not  shown),  we  find  that  the  strong  coupling  greatly 
reduces  the  SNR-improving  effect.  We  may  therefore  draw  another  important  con¬ 
clusion:  there  is  evidently  an  optimal  value  of  k.  If  k  is  too  large,  the  summing 
effect  does  not  work.  If  k  is  too  small,  very  little  signal  will  be  coupled  into  the 
summing  network,  and  the  summed  signal  will  fall  below  the  noise  floor  of  the 
measurement  system. 
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